Quantum Imaging beyond the Diffraction Limit by Optical Centroid Measurements 
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I propose a quantum imaging method that can beat the Rayleigh-Abbe diffraction limit and achieve de Broglie 
resolution without requiring a multi-photon absorber or coincidence detection. Using the same non-classical 
states of light as those for quantum lithography, the proposed method requires only optical intensity measure- 
ments, followed by image post-processing, to produce the same complex quantum interference patterns as those 
in quantum lithography. The method is expected to be experimentally realizable using current technology. 
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PACS numbers: 42.50.St, 42.50.Dv 

It has been suggested that an ensemble of photons, under ap- 
propriate measurements, can be regarded as a single quantum 
object with a much smaller effective wavelength, called the 
photonic de Broglie wavlength [1]. The de Broglie wave- 
length of N photons, each with classical wavelength A, can be 
as small as A /N. It is especially desirable for imaging appli- 
cations to take advantage of the small de Broglie length scale, 
since the resolution of classical optical imaging is limited by 
the size of A, according to the Rayleigh-Abbe diffraction limit 
J2] • A seminal paper by Boto et al. suggests that the -photon 
absorption patterns of certain non-classical iV-photon states 
exhibit de Broglie resolution ||3|]. Subsequent work by various 
researchers on quantum imaging has sought to improve upon 
the original "quantum lithography" technique Jj, |5|, |a [3, |8fl, 
yet all of these proposals still require a multi-photon absorber 
as the detector ||4j, |5[ |g, |8fl or coincidence detection ||7|] . The 
low efficiency of A/-photon absorption and coincidence detec- 
tion, especially for large N, is a significant obstacle to the use 
of quantum technology in real-world imaging applications. 

In this Letter, I propose a quantum imaging method that 
achieves de Broglie resolution without requiring a multi- 
photon absorber or coincidence detection. I shall show how 
the same images as those in quantum lithography can be 
reconstructed simply by intensity measurements, followed 
by image post-processing. Given the availability of high- 
efficiency single-photon detectors, the proposed method has 
the potential to beat the diffraction limit by a large amount 
using current technology. 

A different way of achieving de Broglie resolution is via 
optical-beam displacement measurements, using a split detec- 
tor or homodyne detection, as investigated by Fabre and co- 
workers OJ. The standard quantum limit to the displacement 
uncertainty is on the order of Ax/ y/N, where Ax is the classi- 
cal beam width, while the "Heisenberg" limit is on the order 
of Ax/N 190 ■ The quantum imaging theory I propose here can 
be considered as a generalization of Fabre et al.'s results, and 
therefore presents a unified view of different quantum strate- 
gies that take advantage of the photonic de Broglie wave- 
length. I show that the displacement measurements can not 
only be considered as a Gaussian single-parameter estimation 
technique, but repeated measurements of the parameter can 
also yield complex quantum interference patterns identical to 
those in quantum lithography. The proposed method thus pro- 



vides an accessible way of detecting non-classical and non- 
Gaussian signatures of multi-mode quantum optical states, 
and is expected to become an important tool in the study of 
quantum optics and quantum physics in general. 

I shall first briefly review the configuration-space theory of 
quantum lithography |10y . For simplicity, I consider only free- 
space photons observed at the image plane in one transverse 
dimension, x, in the paraxial regime. Generalization to two 
transverse dimensions and the non-paraxial regime is possi- 
ble HTot [Till , but does not add much insight and the results 
should remain qualitatively the same. Let d(k) and cft(k) be 
the photon annihilation and creation operators, respectively, 
in the transverse-momentum space, with the bosonic commu- 
tation relation [a(k),d^ (k 1 )] — 8(k — k'). First consider a pure 
A/-photon quantum state \N). The multi-photon momentum 
eigenket is \k h ...,k N ) = (N\)- 1 / 2 S t (ki)...d^(k N )\0). The 
momentum wavefunction representation of \N) is then 

<l>{k h ...,k N ) = {ki,...,k N \N)- (1) 

The Rayleigh-Abbe diffraction limit J2l restricts the trans- 
verse momenta of photons to a finite bandwidth: 

2k sin 9 



(fci , . . . , k N ) = for any \k n \> 



(2) 



where sin 9 is the numerical aperture of the optical system and 
sin 9 <C 1 defines the paraxial regime. 

The spatial annihilation operator is defined as A{x) = 
(2n)~ x l 2 j dkd(k) exp(ikx) and can be used to construct 

-photon states with definite positions |jci,...,jcjv) = 
(A r !) _1//2 A t (xi) . . .A^(jcjv)[0). A multi-photon-position posi- 
tive operator-valued measure (POVM) can be defined as 
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with normalization / dx\ . 
tion wavefunction is then 



. dx^Tl(xi , . . . ,Xfj) = 1. The posi- 



\jf(xi,...,XN) = (xi,...,x N \N), 



(4) 



which is the iV-dimensional Fourier transform of 
(j)(ki , . . . Changing the position variables to the 

centroid and relative-position coordinates, defined as 



X ~ N^ x "' 



In 



X, 
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FIG. 1: (Color online). Configuration-space picture of quantum 
imaging for N = 2. The Rayleigh-Abbe diffraction limit restricts the 
momentum wavefunction inside an W-dimensional box \k„\ < ko = 
27Tsin£)//l. The magnitude of the total momentum K = Y^ 1= \^n 
is then limited by \K\ < Nko. The minimum feature size of the 
wavefunction with respect to the centroid coordinate X is then ~ 
l/(Nk )~X/N. 

respectively, a new POVM and a new wavefunction are ob- 
tained: 

?(X,$ 1 ,...,^ 1 )=Ntl(X + $ 1 ,...,X + $ N ), (6) 
f(X,^,...,^ 1 ) = VNw(X + ^,...,X + ^ N ), (7) 

where ^ = — <■;«■ It is not difficult to show that the 
-photon-absorption probability distribution of the -photon 
state p N = \N) (N\ is given by 

(: t N {x) :) = ([P(x)] N [A(x)] N ) « Tr [f (*,0, . . .,0)p N ] 

= |/M,...,0)| 2 , (8) 

which is proportional to the conditional 
centroid probability distribution Pc{x) = 
Tr\T(x,0,...,Q)fa]/ JdxTr\r(x,Q,...,Q)p N ], given that 
all relative positions are zero. The variable conjugate to the 
centroid position X is the total momentum K = Y,n=i The 
bandwidth limit on K becomes \K\ < InNsmQ /X, so the 
minimum feature size of the centroid distribution is on the 
order of A/(2A r sin0). See Fig. [T]for an illustration of the 
configuration-space theory. 

Multi-photon absorption is not the only way of accessing 
the centroid degree of freedom. Consider the intensity cen- 
troid operator, 

X = ^ J dxxA* (x)A (x) , (9) 

which can be measured by spatially resolving intensity mea- 
surements, using, for example, a photon-counting detector ar- 
ray, followed by a calculation of the centroid of the intensity 
pattern, as shown schematically in Fig. [2] The effect of finite 
detector pixel size is to discretize the image plane or equiv- 
alently limit the spatial bandwidth of the photons by virtue 
of the sampling theorem, so the size of each pixel should 
be much smaller than A/(2sin0). If the pixel size is small 
enough, such that the probability of more than one photon 



falling on each is much smaller than one, single-photon de- 
tectors can be used instead. The intensity operator can be 
rewritten as 

X =X1 = X J dx\ . . .dxf]\x\ , . . . ,jqv) (x\ , . . . , jCAf | 
= j dx\ ...dx N (— £x„^ ft (*!,...,*#) 

= JdXd^...d^ 1 Xt(X,^,...,^ l ). (10) 

Thus, measurements of the intensity centroid realizes the 
marginal POVM 

f(X) = f d^...d^ N ^{X,^,...,^ x ), (11) 

and the probability distribution of the intensity centroid mea- 
surements is 

Tr[T(X)p N ]= jd^...d^ 1 \f(X^ u ...^ N _ l )\\ (12) 

which is the marginal centroid probability distribution, here- 
after denoted as p m (X). The minimum feature size of the 
marginal centroid probability distribution is similarly limited 
by the bandwidth of the total momentum < 2nNsin9 / X, 
and can therefore also reach the de Broglie length scale ~ 
X/N. Unlike quantum lithography, which requires an N- 
photon absorber for an A^-fold resolution enhancement, inten- 
sity centroid measurements require only photon counting and 
post-processing to achieve the A^-fold resolution enhancement 
for an arbitrary number of photons. Moreover, as the probabil- 
ity distribution of intensity centroid measurements is marginal 
and not conditioned upon specific values of the relative posi- 
tions, the efficiency of intensity centroid measurements is fun- 
damentally higher than that of multi-photon absorption or any 
other coincidence detection method. This is the central result 
of this Letter. 

In general, the conditional centroid distribution, p c (x), ob- 
tained by multi-photon absorption, and the marginal centroid 
distribution, p m (X), obtained by intensity centroid measure- 
ments, are not the same, except for certain special quantum 
states. One such class of states are the momentum-correlated 
or position-anti-correlated states, of which the "NOON" state 
is a famous example (IH [l(J. These states have small uncer- 
tainties in the relative momenta (k„ = k„—K) and large uncer- 
tainties in the relative positions §„, such that the wavefunction 
can be approximated as 

...,&v-i)«s(*), (13) 

and p c (X) w p m {X) w \g{X)\ 2 . p c (X) and p m (X) become the 
magnitude squared of a complex wavefunction, and can there- 
fore exhibit the same interference patterns as those in classical 
coherent imaging. The momentum-correlated states are also 
the ones that achieve de Broglie resolution, since the band- 
width of the relative momenta K n is negligible, allowing the 
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FIG. 2: (Color online), a) To measure the intensity centroid distri- 
bution, first make an ideal spatially resolving intensity measurement 
(A^ (x)A(x)). This can be done by, for example, an array of photon- 
counting detectors, each with size a -C A/(2sin0). b) The centroid 
X is calculated from the measured intensity pattern, c) The inten- 
sity centroid distribution p m (X) can then be obtained after repeated 
measurements of X, 



bandwidth of the total momentum K to reach its maximum 
value ~ 2?riVsin 9 /X ifioll . Such states can be created, for ex- 
ample, by spontaneous parametric down conversion, and have 
been experimentally demonstrated by D'Angeloef al. ||5|]. The 
time-domain version of momentum-correlated states has also 
been studied theoretically and experimentally by Giovannetti 
and co-workers [12]. 

Another class of quantum states with identical p c (X) and 
p m (X) are the ones with separable wavefunctions as follows: 



/(X,&,...,fr-i)=s(X)A(&,".,&v-i), 



(14) 



of which the quantum Gaussian beams [ 10] and quantum soli- 
tons B13I1 are notable examples. Classical Gaussian beams 
are a special instance of quantum Gaussian beams, and both 
multi-photon absorption and intensity centroid measurements 
of a classical Gaussian beam produce a Gaussian spot with 
a width on the order of X/(VNsin9), as one would expect 
from a classical theory. At the "Heisenberg" limit, the width 
of the centroid distribution of quantum Gaussian beams is on 
the order of X/(Nsm9). One can use the quantum Gaussian 
beam centroid as a "laser pointer" to transmit and reconstruct 
an image in a point-by-point fashion lalialM]- 

Next, I shall generalize the preceding results to quantum 
states with indefinite photon numbers and show that the ef- 
fect of indefinite photon numbers does not significantly affect 



pattern formation by intensity centroid measurements. Define 
a pure state as I 1 ?) = T.n=o^n\N), and denote the wavefunc- 
tion of each Fock-state component by a subscript N. The M- 
photon absorption distribution is 



\Ck 



N=M 

X \\lf N (x,...,X,XM+\ 
M terms 



(Jxm+1 ■ ■ ■ dxtt 
>■■■ ,*n)\ 2 , 



(15) 



which ceases to be the centroid distribution unless N = M. 
This result places stringent requirements on both the gener- 
ation and the detection of non-classical states for quantum 
lithography. 

To properly define an intensity centroid measurement, on 
the other hand, first define the total photon-number operator 
as N = / dxA\x)A(x). Noting that [Nj dxxA\x)A(x)] = 0, 
the intensity centroid operator can be more generally defined 

as 



X= ^ / (/xv.n.viAl.vs 



(16) 



which can be measured by first recording the intensity pat- 
tern, discarding the results when the total photon number is 
zero, and calculating the centroid of the intensity pattern. The 
corresponding POVM is 

f(X)=£ /^i...d§ v _if(X,§i J ...,§v-i) 1 (17) 



T = |0>(0|, Y + JdXY(X) = l. (18) 
The post-selected intensity centroid distribution is then 



Pm(X) • 



N=l 



c N \ 2 dS 1 ...dt N - 1 \f N (x,s 1 ,...,£ 



■N-l, 



(19) 



which is the incoherent superposition of all Fock-state 
marginal centroid distributions. For mixed states, the distri- 
butions are simply the statistical average of pure-state results. 
If a quantum state has a small photon-number uncertainty and 
most Fock components have similar centroid distributions, the 
intensity centroid distribution, being an average of the Fock- 
state centroid distributions, would only be slightly smoothed 
by the incoherent superposition effect, whereas for quantum 
lithography the effect ofN^M can be more drastic |4|, |6J] . 

The centroid measurement can also be used to estimate the 
shift of an optical beam transverse position, due to deflec- 
tion by a mirror in an atomic force microscope, for exam- 
ple. The results concerning the quantum uncertainty of the 
beam displacement obtained by Fabre et al. |@] can be un- 
derstood using the theory proposed here. Making the strong 
mean-field approximation A(x) = A{x) + AA(x), N and X 
can be linearized as N = N + AN and X = X + AX, where 
AN w / dxA* (x) AA(x) + H.c, 



AX ~ I dxxA* (x)AA(x) 



H.c. 



(20) 
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H.c. denotes Hermitian conjugate, and [AN, AX] = 0. AX be- 
comes a quadrature operator that can be measured by homo- 
dyne detection with a local oscillator field <x xA* (x). The the- 
ory presented in this Letter is in fact more general; it shows 
that repeated measurements of X, without the linearization, 
can not only produce a Gaussian spot, but also complex quan- 
tum interference patterns like those in quantum lithography at 
the de Broglie resolution. For position tracking applications, 
such as atomic force microscopy, the proposed method can 
therefore accurately determine the position of a deflecting ob- 
ject even if the position significantly deviates from the mean. 
If the object is quantum, intensity centroid measurements also 
enable the imaging of its position wavefunction at the pho- 
tonic de Broglie resolution. 

Optical loss and imperfect detector efficiency are major is- 
sues for the use of non-classical states for quantum imaging. 
For the specific case of a quantum Gaussian beam, the width 
of p m {X) after the beam propagates through a lossy channel 
can be analytically calculated using the quantum Langevin 
analysis described in Ref. Ill 411 . Assuming (X) = without 
loss of generality, the result is 



(AX 2 ) Z = (AX 2 



Ax 1 



(1-T\e 



(21) 



where (AX 2 ) is the centroid variance, the subscript (z or 0) 
denotes the propagation distance, N z = Noe~ az is the reduced 
photon number, Ax 2 = (N~ l f dxxA' (x)A(x)) is the classical 
beam width squared, given by 



Ax 2 = 



1 



AAk 2 



Rq (l-l/ATp) 2 

JV l-l/(iVb*b) 



(22) 



77 is the detector efficiency, a is the power loss coefficient, and 
Ale ~ 4k sin 9 /X is the root-mean-square momentum band- 
width. Ax and Ak remain constant in a lossy channel UU- 
Rq = 4A^)Afe 2 (AX 2 )o is the initial normalized centroid vari- 
ance, with Ro = 1 at the standard quantum limit (SQL) for 
classical Gaussian beams and Rq = I /Nq at the "Heisenberg" 
limit. In the classical case, Ro = 1, Ax = 1 /(2Ak), (AX 2 ) Z = 
Ax 2 / (t]N z ), and the centroid variance remains at the SQL for 
the detected photon number. Close to the "Heisenberg" limit, 
however, Ax becomes much larger than 1 / (2Ak), and the in- 
crease in centroid variance becomes much more sensitive to 
loss. As long as Rq ^> l/N, quantum imaging by intensity 
centroid measurements is fairly robust to loss. The effect of 
loss for more complex quantum images may be similarly stud- 
ied using the Langevin method, but is beyond the scope of this 
Letter. 



Given the recent progress in the spatial engineering of 
biphoton states IJ, |5[] and photon-counting detector technol- 
ogy, the proposed method is expected to be immediately real- 
izable using current technology, at least for N = 2. For N > 2, 
it is more challenging to create the required non-classical 
states of light. Because of the post-processing requirement, 
the proposed method may not be directly applicable to lithog- 
raphy, while the difficulty of encoding an image in the optical 
centroid degree of freedom limits its practicality in conven- 
tional imaging applications apart from image transmission. 
The proposed method nonetheless greatly simplifies the de- 
tection of multi-mode non-classical optical states, and should 
be useful, at the very least, as a diagnostic tool in the study of 
quantum physics. 
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